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Abstract 

The density and temperature dependence of the nucleonic single particle spectral 
function in symmetric nuclear matter at finite temperatures and densities beyond 
normal nuclear matter density is investigated in a model emphasizing short-range 
correlations and phase space aspects. A simple but self-consistent approach based 
on quantum transport theory is used. In particular we consider the density and 
temperature regime occuring e.g. during the core collapse in a supernova explosion 
and the subsequent formation of a neutron star. Mean-field effects are incorporated 
by a Skyrme type potential. 
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1 Introduction 



An interesting but only partially explored problem of nuclear many-body the- 
ory are contributions from short-range correlations to the binding properties 
of nuclear matter and finite nuclei. Of particular interest for various regions 
of nuclear physics ranging from heavy ion collisions to the physics of neutron 
star formation is the density and temperature dependence of short-range cor- 
relations in nuclear matter. In heavy ion collisions like those planned at the 
Compressed Baryonic Matter (CBM) facility at GSI densities of several times 
the nuclear matter saturation density po = 0.16/ fm 3 are reached. Such condi- 
tions are also encountered in astrophysical scenarios. When a neutron star is 
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formed in a supernova explosion temperatures of several tens of MeV and den- 
sities far beyond normal nuclear matter density arise. The neutrinos emitted 
from the nascent neutron star are an important signal that gives information 
about the ongoing cooling processes. Nucleon correlations, however, affect di- 
rectly the neutrino opacity in dense nuclear matter. The neutrino cross sections 
can be significantly reduced due to the presence of correlations in the medium 
[1,2]. Therefore it is important to take nucleon correlations into account when 
interpreting the neutrino spectra. 

So far there has been no consistent investigation of the properties of the nu- 
cleon spectral function that covers the full temperature and density range 
of those scenarios. In this work we extend the approach to short-range cor- 
relations in nuclear matter presented initially in [3,4] to the more general 
cases of arbitrary temperatures and densities. Using direct relations between 
correlation functions and collision integrals a simple but fully self-consistent 
model for calculating the nucleon spectral function in an iterative procedure 
was constructed in refs. [3,4]. The former calculations, however, have been 
restricted to zero temperature and normal nuclear matter density. The pre- 
viously obtained results are in surprisingly good agreement with much more 
sophisticated and numerically extremely involved, "state-of-the-art" calcula- 
tions utilizing the full machinery of modern nuclear many-body theory [5,6]. 
In particular, the significant population of high-momentum states in the nu- 
cleon momentum distribution - which is interpreted as an overall measure for 
short-range correlations - is described rather well. 

In this paper, we reconsider this model and extend it to finite temperatures 
and densities beyond the nuclear matter equilibrium point. Different from [3,4] 
this requires a reliable description of the static mean-field self-energy which 
in the former work could most conveniently be absorbed into a redefinition of 
the chemical potential. In the present case an explicit treatment of the density 
and momentum dependence of the mean-field self-energy is necessary since the 
chemical potential is temperature and density dependent. Calculations for the 
equation of state for dense matter and finite temperatures have been made 
using microscopic Hamiltonians, see e.g. [7,8]. In our model we incorporate 
the static properties of the mean-field by an empirical, Skyrme-type energy 
density functional [9,10,11]. By this choice thermodynamical consistency is 
also ensured by fulfilling the Hugenholtz-van Hove theorem [12] at mean-field 
level. This approach gives rise to a potential that is independent of energy but 
depends on momentum, temperature, and density. However, because of the 
Skyrme-typical restriction to terms of at most quadratic order in the single 
particle momenta the momentum dependence can be re-expressed immediately 
by an effective mass in the kinetic energy term. 

Calculations similar to ours have been done by Alvarez-Ruso et al. [13] for 
densities from half to twice p and temperatures up to 20 MeV. In their "semi- 
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phenomenological" approach they have evaluated second order diagrams tak- 
ing the magnitude of the NN interaction from experiment. However, they did 
not perform their calculations self-consistently and did not consider long-range 
contributions from the mean-field. Note that the basis of their work was also a 
model for cold nuclear matter [14]. Self-consistent approaches to the nucleon 
spectral function in nuclear matter have been studied theoretically before, 
e.g. by Aim et al. [15] and in ref. [16]. More recently, calculations using a 
self-consistent T-matrix approach have been presented by Bozek [17,18]. In 
[17] the temperature dependence of the spectral function was investigated up 
to T = 20 MeV, but the density dependence has so far not been investigated 
by Bozek. 

In this work we go far beyond the temperature and density ranges of [13]- 
[18]. To study the role of short-range correlations in supernovae and heavy 
ion collisions at CBM we consider temperatures up to 70 MeV and densities 
from po to 3po- in section 2 we briefly outline our model. The results of the 
calculations are presented in section 3. They are compared to other models (at 
T = and finite temperatures) and we discuss the temperature and density 
dependence of the short-range correlations in nuclear matter. A summary is 
given in section 4. 



2 The Model 



The model we present here was introduced in [3,4] for cold (T = 0) nuclear 
matter at saturation density po- The underlying Green's function formalism 
is described in much detail in [19,20]. Thus we restrict ourselves here to a 
short summary of the concept. We show in more detail how to modify the 
original model [3,4] for the application to arbitrary densities and temperatures. 
Throughout this paper we consider only symmetric nuclear matter. 

On the two-particle-one-hole (2plh) and one-particle-two-hole (lp2h) level, 
respectively, the collisional or polarization self-energies S > and S < are given 
by: 



> r d 3 p 2 duj 2 d 3 p 3 du 3 d 3 p 4 duj A 4 4 \TTl2i 

s = 9 J 727F7^F7^T ( } {p + P2 " n " Va) lM 1 

xg^(uj 2 , p 2 )g^(u 3 , p 3 )g^(uj 4 , p 4 ), (1) 

where g = 4 is the nucleonic spin-isospin degeneracy factor and \Ai\ 2 is the 
square of the nucleon-nucleon scattering amplitude averaged over spin and 
isospin of the incoming and summed over spin and isospin of the outgoing 
nucleons. We leave the structure of M. open for the moment and come back 
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to it later. The correlation functions g > and g K in eq.(l) can be rewritten in 
terms of the spectral function a(u,p), 

g < (uj,p)=ia(uj,p)f(uj,p), 

g > (u,p) = -ia(u;,p)(l - f(u,p)), (2) 

with the phase space distribution function f(u,p). In thermal equilibrium / 
is just the Fermi distribution function depending on the energy and, through 
the chemical potential u>f, on the density of the system, 



/M = 1 + e (I^)/T - ( 3 ) 

The previous calculations were greatly simplified by the fact that at zero 
temperature the Fermi distributions reduce to step functions allowing for con- 
siderable simplifications of eq.(l) by taking advantage of closed form analytic 
expressions [3,4]. 

The single particle spectral function in (2) is explicitly given as 



aiuj.p) — 5 , 4 

(w - £ - £ mf - ReS«*(a;,p)) 2 + ip( w> p)' 

where ReS ret is the real part of the collisional self-energy and T is the width 
of the spectral function. T is directly related to the imaginary part of the 
retarded self-energy, 



T(u,p) = -2lm^ ct (uj,p) = i{T>{u,p) - E<(u,p)). (5) 

Note, that eq.(4) also includes the mean-field self-energy S mf = S mf (p, p, T) 
which is not contained in E^. The purely hermitian S mf contains the contribu- 
tions to the self-energy arising from ground state tadpole diagrams describing 
mainly long-range correlations. 

We can rewrite the first term in the denominator of eq.(4) by expanding 
S mf (p, p, T) w S mf (k F , p, T) + (p 2 - 4)E mf '(A; F , p,T) + ... up to 0(p 2 ). The 
momentum dependent part of S mf can then be absorbed into the kinetic en- 
ergy term by introducing a density and temperature dependent effective mass 
= — (1 + 2mS mf '), as successfully used, e.g., in Skyrme density functionals 
[9,10,11]. We get: 



2 2 

|- + S mf (p,p) + ReS rct (^,p) = + U cS (p) + Re£ r( >,p), (6) 
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where the effective potential U e g is the momentum independent part of E . 
We will discuss the role of U e s and ReE ret in more detail below. 

The equations (l)-(5) describe a self-consistency problem. Typically, such 
problems cannot be solved directly but iterative procedures are the appropri- 
ate methods of solution. Hence, the spectral function is obtained iteratively 
by the following procedure: Starting with an initial guess for a(u,p) in the 
integrals of eq.(l) an approximation for the width T is found. Inserting this 
T(uj,p) into eq.(4) yields an improved expression for the spectral function that 
is inserted into (1) again. This procedure can be repeated until convergence 
is achieved. In the language of Feynman diagrams the solution of this self- 
consistency problem corresponds to the non-perturbative summation over a 
whole class of diagrams by reinserting the full in-medium Green's functions 
into the self-energies. 

As we have seen above there are two real contributions to the self-energy. One 
is the mean-field self-energy that is purely real and does not affect the width 
directly. The other contribution is the real part of the collisional self-energy 
ReE rct . It is related to the imaginary part of the polarization self-energy (i.e., 
the width of the single particle spectral function) by a dispersion relation, 
e.g. discussed in [4]. There, the importance of the dispersive contributions was 
investigated by performing calculations with and without ReE rct . Although 
analyticity is not conserved when ReE ret is neglected [3] it was found that 
overall the results on spectral functions did not change very much. Taking 
advantage of this finding we will ignore the dispersive part for the present 
investigations, thus simplifying the calculations considerably. For a full treat- 
ment of ReE ret we refer to [3,4]. 

More important for our present purpose of studying the density and tem- 
perature dependence of short-range correlations are the mean-field potential 
U c fi = U e g(p,T) and the effective mass m* = m*(p). The net effect of the 
potential U e g is a temperature and density dependent energy shift of the pole 
structure of the single particle propagators and hence the spectral functions. 
For the investigations in [3,4] at T = and p = p such a shift could be ab- 
sorbed into a re-definition of the chemical potential. This is no longer possible 
when we look for the density and temperature dependence of our results to 
which long-range correlations contribute. 

The temperature and density dependence of the single-particle potential in 
nuclear matter has been studied before microscopically, ranging from self- 
consistent Brueckner and Brueckner Hartree-Fock approaches e.g. in [15,16] 
to the variational calculations with 2- and 3-body forces by the Urbana group 
[7,8]. In a different context, the density and momentum dependence of the 
variational mean-field was studied in [7] . A re-analysis of these rather involved 
mean-field results by a second order polynomial in k 2 — k 2 F leads to an almost 
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perfect description of the momentum dependence with an effective mass of 
m* jm = 0.69. This confirms the effective mass approach utilized in eq.(6). 
These findings also confirm the use of a schematic description of the mean- 
field part, in particular so because the static parts of the self-energies are 
cancelled to a large extent in dynamical quantities like propagators. Hence, 
we decide to describe the mean-field in a Skyrme model [9,10,11] for the sake 
of a most transparent and simple approach. 

From the Skyrme energy density functional in standard notation [9,10,11] U c s 
and m* in infinite symmetric nuclear matter are obtained by 



U eS (p) = h oP + ^(3tx + 5t 2 )r + ^t 3 p 2 , (7) 

i = i + i> i + 5t ^ (8) 

where to, h, £2, and £3 are parameters of the Skyrme interaction [9,10,11]. The 
parameter set we have used here is given in Table 1. r is the (quasi-particle) 
kinetic energy density, 



7T 



dp 



p 

; (p 2 -fcf.)/2 m *T + [' 



(9) 



with kp = 2m*(uJ F - U eS ). 

It is well known that nuclear matter ceases to be bound beyond a critical 
density p c {T) where the repulsive kinetic energy term overrules the nuclear 
attraction. For processes inside a supernova this is, however, not a crucial 
problem. The core collapse is not driven by nuclear, but by gravitational forces. 
Eventually, the repulsive nuclear potential even stops the collapse when several 
times nuclear matter density is reached. 

For a complete definition of the approach we also have to specify the scattering 
amplitude M. in the collision integrals (1). Taking into account the observation 
of Lehr et al. [3,4] that the short-range correlations primarily depend on global 
properties of the interactions on ranges of the (repulsive) vector meson masses 
and the available phase space we will also use a constant matrix element in 
our calculations, corresponding to an effective zero-range contact interaction. 

The value of M. will be chosen independent of temperature and density of 
the system. It is noteworthy that until now a description of short-range cor- 
relations including a fully self-consistent description of the interactions is still 
pending, although first attempts in that direction have been made. In such 
calculations, e.g. by the Rostock group [15] or, more recently, by Bozek [17,18] 
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significant variations of the interactions with temperature are not obtained. 
Accepting that short-range correlations are due to short-range interactions, a 
strong dependence of the interaction matrix element is not very likely, either. 

In order to enforce the convergence of integrations involving Ai and extending 
over the space of 2plh and lp2h configurations, respectively, we introduce a 
form factor such that the collision integrals in eq.(l) are cut off at large u. 
Lehr et al. had introduced a form factor that acts on the total energies and 
momenta of the collisions corresponding to eq.(l). This form factor did only 
affect the results at uo > uo F . For the occupied states below uj f the results 
remained unchanged. Due to the density dependence of the chemical potential, 
however, this form factor is not a sensible choice for our calculations. 

Different to the approach in [3,4] we use here a form factor acting on the t- 
channel and suppressing processes with large energy and momentum transfers 
(see [21] for details): 



where u t = uj — uo^ and Pt = p — P3- The mass m M ~ 0(m p ) corresponds to 
the typical mass of an exchanged meson in nucleon scattering processes and 
A defines the scale of the cutoff. Note that F has no pole on the real axis. 



3 Numerical Details and Results 

We have calculated the nucleon spectral function in nuclear matter at temper- 
atures from 10 to 70 MeV and at densities p — p • • • 3p (Po = 0.16 fm -3 ). In 
the calculations we used the Skyrme force type SIII' to generate a mean-field 
potential and an effective mass m* (cf. [9,10,11]). The corresponding set of 
parameters is given in Table 1. In infinite nuclear matter at zero temperature 
one finds for this choice of parameters a binding energy of —16.1 MeV and a 
mass ratio of m* jm = 0.76 at a Fermi momentum of k F = 1.29 fm -1 (p = p ). 
For the mean-field potential a value of U e ff = —61.6 MeV is found. Compared 
to the variational results of Wiringa [7] the saturation properties are slightly 
different with an equilibrium point shifted to a somewhat lower density. The 
larger effective mass indicates a weaker momentum dependence of the SIII' 
interaction. 

The averaged scattering amplitude was adjusted to the value 1 ( |_Axt | 2 ) x / 2 = 
1 This value differs from the one used in [3,4] by the factor m/m* « 1.4. 




(10) 
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309MeVfm 3 so that the results of the many-body calculations of Benhar et 
al. for T = [5] are reproduced. The cutoff parameter A of the form factor 
(10) was set to 1.2 GeV and the mass m M to 600 MeV. To verify that this is a 
reasonable choice for the parameters we compare the on-shell width r on (c<j) = 
T(uj,p(uj)), where p(uo) = (2m*(w-[/ cf f)) 1 / 2 , from our calculation to the results 
of Benhar et al. [5] and Baldo et al. [6]. In Fig. 1 the on-shell widths are shown 
for normal nuclear matter density and T = 0. As it can be seen, our result 
is in good agreement with the sophisticated many-body calculation of Benhar 
et al. and supports the choice of the parameters. Note that the discrepancies 
between the results of Baldo and Benhar are explained in [5] by the different 
NN interactions that are employed. 

In the following we will show the results of our calculations and investigate 
the density and temperature dependence of the short-range correlations. First, 
we present cuts of the nucleon spectral function at a constant momentum 
to illustrate the general role of temperature and density. Fig. 2 shows the 
spectral function for temperatures of 0, 10 and 50 MeV at normal nuclear 
matter density. The effect of higher temperatures is clearly visible. The gap at 
the Fermi energy gets smeared out for T = 10 MeV and has almost vanished 
at T = 50 MeV. This phenomenon is directly related to the fact that for 
non-vanishing T - in contrast to a system at zero temperature - states at 
energies above up are populated. The states at the Fermi surface cease to be 
quasi-stable. 

In Fig. 3 the spectral function is shown for three different densities at a con- 
stant temperature of 10 MeV. Since the chemical potential depends on the 
density the three curves are shifted along the energy axis. Note that the rela- 
tion fj, ~ p 1//3 for the free Fermi gas does not hold here as it can be directly 
seen for T = 0, where \i = r + Re(E(/ci?)). One should not get mislead by 
the observation that the on-shell peaks at lower densities are broader than 
the peak at the highest density. This effect is due to the fact that the on-shell 
energy lies closer to the Fermi energy for the highest density. This induces a 
stronger suppression of the on-shell peak than in the other cases. 

Figs. 2 and 3 give a first impression of the spectral functions. To investigate the 
temperature and density dependence of the short-range correlations, however, 
we use the better suited on-shell width of the spectral function, T on (uj). In 
Fig. 4 we show our results for T on (uj) at temperatures of 10 and 20 MeV and 
densities of p and 2 p . Also shown are the results of Alvarez-Ruso et al. 
[13]. As we have already seen in Figs. 2 and 3 the width of the states at 
the Fermi surface does not drop to zero for finite temperatures. Increasing 
the temperature leads to larger values of r on in the vicinity of the Fermi 
surface while the size and the shape of r on do not change significantly at 
higher energies. Increasing the density shows the opposite effect. r on is left 
unchanged in the vicinity of the Fermi surface while the slope at high energies 
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increases. Compared to [13] our results seem to be shifted by approximately 
5 MeV to higher values of r on . The shift can be explained by the differences 
in the NN interactions that were used in both calculations. The shape of the 
curves is in good agreement for u — u F < 150 MeV. At higher energies the 
form factor that was used in our calculations affects the shape of r on , e.g. 
leading to the observed curvature in our results. 

Qualitatively our results (for p = p Q and T < 20 MeV) are similar to those 
obtained in the self-consistent T-matrix approach by Bozek [17,18] but on a 
quantitative level significant differences are found. However, the deviations 
range within the uncertainties introduced using different nucleon-nucleon in- 
teractions in fully microscopic calculations. This is evident by comparing the 
results in [18], based on the Paris NN potential, e.g. to the variational results 
of Benhar et al. [5], obtained with the Argonne NN potential. The reason for 
these deviations is found in the fact that by the way of deriving these empirical 
interactions from NN data they are guaranteed to agree in their predictions 
of on-shell quantities, as e.g. free space NN scattering phase shifts, but might 
differ in the off-shell region. While we work with a constant coupling in com- 
bination with a t-channel form factor Bozek tends to use strongly momentum 
dependent interactions like the Paris potential (see [18] for details). Similar 
observations are discussed in [18] when results of different approaches - in- 
cluding [3] and [17] - are compared. Note that the same effect has already 
shown up in Fig. 1 for the results of Benhar [5] and Baldo [6] . 

The temperature dependence of the on-shell width at various densities is illus- 
trated in Fig. 5. The effect of filling the gap at the Fermi surface with increasing 
temperature is clearly visible. At all densities one observes a levelling off in 
the energy dependence of T oi1 (uj) at the highest temperatures indicating the 
gradual changes in the available phase space in a heated system. 

The behavior seen in Fig. 5 deserves a closer inspection. Before doing so we 
should take into account that the chemical potential uof depends strongly on 
the nucleon density, as was shown in Fig. 3. Thus one has to be careful at 
which energy u and momentum p(u), respectively, the widths for different 
temperatures and densities are compared. This problem is avoided by aver- 
aging the on-shell width over the occupied states and compare this (r on ) for 
the different p and T. Fig. 6 shows the overall behavior of (r on ) for the full 
temperature and density range that was covered in our calculations. Note that 
one gets a very similar picture when T on ((uj)), i.e. the on-shell width at the av- 
eraged energy of the occupied states, is displayed instead. The general features 
of this figure are thus universal. It can be seen in Fig. 6 that the on-shell width 
ranges from 10 MeV at low temperatures up to 40 MeV at high temperatures. 
(r on ) rises almost linearly with the temperature at medium and high densities. 
At normal nuclear matter density the dependence is only linear for tempera- 
tures below Ki 45 MeV and gets weaker at higher temperatures. The density 
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dependence of (r on ) that can be derived from Fig. 6 has a more complicated 
structure. As long as the density is low (p < 2p ) (Ton) does increase with the 
density just as one would expect naively. At densities above 2p , however, the 
density dependence becomes very weak - in particular for temperatures above 
30 - 40 MeV. Since we take the on-shell width as a measure for short-range 
correlations this means that the correlations saturate at densities above two 
times normal nuclear matter density. 

The explanation for the saturation is rather simple. Pauli-blocking suppresses 
nucleon collisions (the source of the correlations) with final states below the 
Fermi surface. Since the chemical potential depends strongly on the nucleon 
density more (kinematically allowed) final states are blocked at higher p. At 
high p there are also more collision partners (initial states) for a particle with 
given energy and momentum but this cannot compensate the suppression from 
Pauli-blocking. To understand this one has to look at the total collision rate 
for a nucleon with given u and p with the nucleons of the medium (cf. eq.(l)). 
There are two final states that are affected by Pauli-blocking but there is only 
one initial state from below the Fermi surface. Thus the effect of increased 
Pauli-blocking is more important than the gain of initial states. 

To verify that this explanation is correct we have made additional calculations 
in which Pauli-blocking was switched off. The result is shown in Fig. 7. It can 
be clearly seen that without Pauli-blocking there is a linear dependence of 
(r on ) on both, temperature and density, over the full range covered by our 
calculations. This is exactly the picture one would naively expect when Pauli- 
blocking is not taken into account. Thus Pauli-blocking is in fact the correct 
explanation for the saturation of the short-range correlations at high densities. 



4 Summary and Conclusions 

In this paper we have presented an approach to short-range correlations in 
nuclear matter at finite temperatures and high densities. Based on the work 
of Lehr et al. [3,4] we have constructed a simple but self-consistent model that 
goes beyond the quasi-particle approximation. The interactions between the 
nucleons were described by a constant matrix element in combination with a 
form factor that limits the energy and momentum transfer in the collisions. In 
[3,4] the averaged coupling constant has led to excellent results for a system 
at zero temperature. We have not calculated the mean-field contributions self- 
consistently within the model. Instead, a Skyrme-type interaction model was 
incorporated to assure a realistic thermodynamical behavior of the chemical 
potential. 

We have calculated the spectral function of nucleons in symmetric nuclear 



10 



matter at temperatures from 10 to 70MeV and densities from po to 3po- This 
has been the first time that this large regime is investigated consistently within 
the same model. It was shown that our results are in good agreement with 
other calculations at temperatures up to 20 MeV and densities up to 2 p . The 
width of the spectral function was used to get informations about the temper- 
ature and energy dependence of the short-range correlations. We have found 
that the correlations scale approximately linear with temperature. At high 
densities, however, the correlations saturate, an effect that can be explained 
by Pauli-blocking. 

The model we have presented here is open for further improvements. Using sep- 
arate spectral functions for protons and neutrons and introducing an isospin 
dependent coupling it would also be possible to investigate isospin asymmetric 
nuclear matter. 
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to 

(MeV • fm 3 ) 


h 

(MeV • fm 5 ) 


t 2 

(MeV • fm 5 ) 


h 

(MeV • fm 6 ) 


x 


W 
(MeV • fm 5 ) 


-1133.4 


395 


-95 


14000 


0.49 


120 



Table 1 



The parameters of the Skyrme force SIII' [9] used in the calculations. For complete- 
ness this table also contains the spin-interaction parameter xq and the spin-orbit 
parameter Wq. 
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Fig. 1. The on-shell width r on (u;) as a function of uj — ojp for T = and normal 
nuclear matter density. The solid line shows the result of our calculation, the dashed 
line the result of Benhar et al. [5], and the dotted line the result of Baldo et al. [6]. 
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Fig. 2. The nucleon spectral function at T = MeV (solid line), T = 10 MeV (dashed 
line) and T = 50 MeV (dotted line) for normal nuclear matter density. All cuts were 
made at a constant momentum of p = 0.45 GeV. 
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Fig. 3. The nucleon spectral function for normal nuclear matter density (solid line), 
a two times higher density (dashed line), and a three times higher density (dotted 
line) at T = lOMeV. All cuts were made at a constant momentum of p = 0.45 GeV. 
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Fig. 4. The on-shell width T ou (uj) as a function of lo — <jjp for two temperatures and 
densities. The solid lines show the results of our calculations, the dashed lines show 
the results of Alvarez-Ruso et al. [13]. 
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Fig. 5. The on-shell width r on (u;) as a function of uj — u>f for different temperatures 
and densities. Results for T = 10, 30, 50, 70 MeV are indicated by solid, dashed, 
dotted and dashed-dotted lines, respectively. 
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Fig. 6. The on-shell width (r on ) in the (T, p)-plane (T in MeV). For each temper- 
ature and density r on has been averaged over the occupied states. 
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Fig. 7. The on-shell width (r on ) in the (T, p)-plane (r in MeV) when Pauli-blo eking 
is switched off in the calculations. For each temperature and density r on has been 
averaged over the occupied states. 
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